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ABSTRACT 

Aims. We report the results of a survey of 442 planetary nebulae at 30 GHz. The purpose of the survey is to develop a list of planetary 
nebulae as calibration sources which could be used for high frequency calibration in future. For 41 PNe with sufficient data, we test 
the emission mechanisms in order to evaluate whether or not spinning dust plays an important role in their spectra at 30 GHz. 
Methods. The 30-GHz data were obtained with a twin-beam differencing radiometer, OCRA-p, which is in operation on the Torun 
32-m telescope. Sources were scanned both in right ascension and declination. We estimated flux densities at 30 GHz using a free-free 
emission model and compared it with our data. 

Results. The primary result is a catalogue containing the flux densities of 93 planetary nebulae at 30 GHz. Sources with sufficient 
data were compared with a spectral model of free-free emission. The model shows that free-free emission can generally explain the 
observed flux densities at 30 GHz thus no other emission mechanism is needed to account for the high frequency spectra. 

Key words. ISM: Planetary Nebulae: general - Radio continuum: general 
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1. Introduction 

The planetary nebula (PN) phase in the evolution of low mass 
stars lasts only about 10 4 years. It begins once the central star 
reaches an effective temperature of 20 000 K and ionises the shell 
of material developed during Asymptotic Giant Branch (AGB) 
evolution. The end of this phase is defined by termination of 
nuclear burning in a thin outer shell of the star, and then rapid 
dispersal of the nebula. At radio wavelengths, planetary nebulae 
emit continuum radiation in a free-free process and are among 
the brightest Galactic radio sources. Their radio flux densities 
do not suffer the high levels of extinction present in the opti- 
cal regime. Since most planetary nebulae occupy Galactic plane 
where extinction is high, radio detections and radio flux density 
measurements are important. 

Planetary nebulae are mostly compact sources because they 
are distant or intrinsically small. Their relatively strong and sta- 
ble radio emission makes them good candidates for calibration 
sources. 

Many radio continuum observations of planetary nebulae 
have been made at 1.4, 5 and 14.7 GHz (seelCon don & Kaplan 
1991 lAcker et alJ[l99i IXaauist & Kwokll 19901: iMilne & Allei 
19821) . There is still only limited data at frequencies above 
30 GHz and most of those observations were obtained with in- 
terferometers such as the VLA; relatively little comes from sin- 
gle dishes. To extend the spectral range and to make total flux 
density measurements, we have used the new One Centimetre 
Radio Array prototype (OCRA-p) receiver to observe planetary 
nebulae at 30 GHz. This receiver is mounted on Toruri's 32-m ra- 



dio telescope and is described in section [2] OCRA-p receiver is 
outlined in detail bv lLowd (120051) . The survey for planetary neb- 
ulae was one of the first successful observations made using this 
system (together with measurements of flat-spectrum sources by 
lLowe et al. 20 071 and observa tions of the Sunayaev Zel'dovich 
effect by Lancaster et alj|2007l) . The purpose was to make a list 
of high frequency calibrators, which can be used to support sky 
surveys and to test the emission mechanisms in order to evalu- 
ate whether or not spinning dust plays an important role in PN 
spectra. 

Our new survey of planetary nebulae brought detections of 
93 sources at 30 GHz out of 442, of which selection criteria are 
described in section|3] The observing techniques and data reduc- 
tion process are described in section|3] The results are described 
in section [4] The comparison of flux densities with the fre e-free 
emission model proposed by Siodmiak & Tvlenda (2001) is de- 
scribed in section[5] Section|6]contains the final conclusions. The 
measured 30-GHz flux densities of all detected PNe are given in 
Tabled 



2. OCRA-p 
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OCRA -p is a 2-element prototype for OCRA dBrowne et al.l 
2000) - a 100-element array receiver. The OCRA-p receiver was 
constructed at the University of Manchester and was funded by 
a EU Faraday FP6 grant together with the Royal Society Paul 
Instrument fund. It has been mounted on the Torun 32-m ra- 
dio telescope owned by the Centre for Astronomy, Nicolaus 
Copernicus University, Torun in Poland. 

The basic OCRA design was based on the prototype 
demonstrator for the Planck Low Frequency Instrument 
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Tabl e 1. OCR A-p nominal system specification (from 
lLowe et al.ll2007h 



Frequency range 
System temperature 
Individual beam FWHM 
Beam spacing 



27 - 33 GHz 
40 K 

1.2 arcmin 
3.1 arcmin 



dMandolesi et alJ 120 00) and is similar to the K-band rec eivers 
mounted on the WMAP spacecraft dJarosik et all 120031) . The 
nominal system specification is presented in Table 1 . OCRA-p 
has two closely spaced feeds in the secondary focus of the 32 m 
telescope and thus there is very similar atmospheric emission 
in each beam. The basic observing mode with this radiometer 
involve switching continuously between the two horns and be- 
tween two states of a 277 —Hz phase switch (0 and 180°), located 
in one arm of the receiver. The switching frequency is much 
greater than the characteristic frequencies of atmospheric and 
gain fluctuations and those effects are dramatically reduced by 
taking differences between switch states and channels. The re- 
sult, so-called double difference, is the observable used for data 
reduction. Mor e details and an ill ustration of our scan strategy 
can be found in lLowe et alJ (120071) . 

The full 100-beam OCRA receiver could be ideal for study- 
ing radio sources at higher frequencies, especially it is ideally 
suited to fast sky surveys. 

3. Observations and data reduction 

The observations were made over a period of 1.5 years, starting 
in December 2005. For a survey we selected 442 planetary neb- 
ulae from the Strasbourg-E SO Catalogue of Galactic Planetary 
Nebulae (Ack er et al Jl 19941) having 5 GHz observations and de- 
clinations above -15". 

A single flux density measurement consisted of two orthogo- 
nal scans with the first one made in elevation. Real-time software 
fitted a Gaussian profile to the observed beam shape and the mea- 
sured position offset was used to update the telescope pointing 
corrections for the subsequent azimuthal scan. These corrections 
of pointing accuracy led to reliable estimate of flux density using 
the azimuthal scans. The first step of the data reduction was to fit 
a quadratic background baseline, subtract it from the data and fit 
double Gaussians to the azimuthal scan to ensure that the source 
is passing through each beam in turn. The quality of the data 
was strongly dependent on the weather conditions. During bad 
weather large asymmetries would be observed in the fitted am- 
plitudes of the two beams. Data with differences in amplitudes 
exceeding 10% were not used. To fix the flux scale, NGC 7027 
was used as the calibrator for the whole survey. It is the brightest 
planetary nebula in the radio sky due to its proximity and youth . 
It is a compact nebula of 10x13.5 arcsec ( Athe rton et alj[T979h 
observed regularly at the VLA for 25 years and can be taken to 
be unresolved when observed with the 1 .2 arcmin beam of the 
32 m telescope. Those features make NGC 7027 a good prime 
calibrator over a wide frequency range. Fortunately, at Toruh's 
latitude the source is circumpolar so it could be observed at any 
time during the survey. 

On human time scales the flux evolution of NGC 7027 is 
slow but it is measurable. iPerlev et"ai] d2006l) suggest that the 
observed smooth decline of PNe flux densities at high frequen- 
cies is due to the reduction in the number of ionising photons 
with time. This is caused by evolution of the central star, whose 
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Fig. 1. Histogram of the distribution of measured planetary neb- 
ulae flux densities at 30 GHz. 



temperature is inc reasing and luminos ity is slowly decreasing. 
Measurements bv lZiilstra et al.l (|2008) show that the flux den- 
sity of NGC 7027 is decreasing at a rate of -7.79mJy/yr at 
30 GHz. Using this secular decrease in flux density, along with 
a flux density of 5433.8 mjy at 2000.0 based on VLA measure- 
ments in seven bands, we calculate a flux density for 2006.0 of 
5.39+0.28 Jy. The data for all observed nebulae were scaled to 
this value. 

The final results were corrected for gain-elevation curve of 
the telescope with elevation (a factor of up to 20%, depending 
of elevation). The atmospheric absorption was also taken into 
account. The final flux density uncertainty arising from double 
Gaussian fitting, calibrator uncertainty and all other corrections 
is about 8%. 



4. Results 

The final results are presented in Table [2] The first column 
contains the SESO IAU Planetary Nebula in Galactic coordi- 
nates (PNG) name Ul.l+bb.b, where 111.1 is Galactic longitude 
and bb.b is Galactic latitude. The second column contains 
the normal descriptive designation. Coordinates for J2000, 
taken from I Acker et ail d 19941) . are given in columns 3 and 
4. Columns 5 and 6 contain the measured flux densities at 
30 GHz and calculated uncertainties. The values marked by * 
are sources with optical angular size larger that 1/4 OCRA-p 
beam, while those marked by A are corrected using adopted 
angular diameter which can be found only for some sources 
(see section |5j. Table [2] contains data for 93 sources. Figure 7 
presents PN modeled radio spectra (see section|5]l together with 
flux densities at 1.4, 5, 14.7 GHz. 

The distribution of detected flux densities in presented in 
Figure Q] It is seen that there is an instrumental limit on the 
detectability of sources around 20 mjy. The exact value depends 
on the weather conditions at the time of observations. 



5. Analysis of flux density results 

The first result of the survey is a new catalogue of 93 Planetary 
Nebulae flux densities at a frequency of 30 GHz. We have plotted 
in Figure [2] a histogram for the distances of all PNe surveyed in 
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Fig. 2. Histogram of the distribution of planetary nebulae dis- 
tances. The plot presents distances to all planetary nebulae ex- 
amined in the survey and the ones detected by OCRA-p are high- 
lighted. 



this work, and marking those detected by OCRA-p. The data was 
taken from the Strasbourg-ESO Catalogue of Galactic Planetary 
Nebulae. Histogram shows that we were able to measure plan- 
etary nebulae flux densities at distances out to about 5 kpc and 
that the distributions of detected and non-detected PNe in our 
data are similar. 

Planetary nebulae are strong infrared sources and it is 
known that they occupy cert ain characteristic po sitions on the 
IRAS colour-colour diagram. iPottasch et all (fl988l) made IRAS 
measurements of over 300 planetary nebulae in 4 bands: 12 jjm, 
25/im, 60 fim and lOOyum. We used the first three to construct 
a useful colour-colour diagram. Figure [3] shows this diagram 
for all examined sources and the ones detected with OCRA-p. 
It can be seen that planetary nebulae measured by OCRA are 
concentrated close to the black-body (120-180K) line in the 
plot. This is due to the evolution of planetary nebula. As the dust 
cools with time, the dust colours will redden whilst the nebulae 
becomes less opaque in the radio region, resulting in an increase 
in radio emission. Thus observed radio emission provides useful 
hints on planetary nebulae evolution. 

Planetary nebula emission at 30 GHz is thought to be due to 
free-free emission. We want to check whether this mechanism 
can fully explain the measured flux densities since there have 
been claims that the excess emis sion has been seen in plane- 
tary nebulae (Casassus et al. 2004) . A model of planetary nebul a 
free-free emission was taken from lSiodmiak & Tvlendal (12001). 
These authors analysed the flux densities of 264 planetary neb- 
ulae with reliable measurements at 1 .4 and 5 GHz and known 
diameters (with radio diameter defined by radio contours at 10% 
of the peak flux density, or the mean value of optical and radio 
diameter). Four simple models were proposed and examined. 
The most reasonable approach is to consider a model consist- 
ing of two regions with different optical depth. Each source is 
then modelled a dense region which fills a solid angle £Q char- 
acterised by optical depth t v , and a less dense region filling up 
(1 -£)£2 described by optical depth er v . With such an assumption 
one can model the spectra of all PNe, including the small frac- 
tion of the population which are homogeneous and symmetric. 



The observed flux density is then expressed as: 
2v 2 kT 

= —r^ ((1 - e^)&l + (1 - e~^) (1 - £)Q) , 



(1) 



where the optical thic kness depends on the frequency v roughly 
as (see lPottaschlll984l) : 



Tv = To(v/v ) 



-2.1 



(2) 



from 
from 



to and vo (=5 GHz) are the reference optical thickness and fre- 
quency. Optical thickness at 5 GHz is found by solving the first 

expression with known flux densit y at this frequency. 

The best fit values found by [Siodmiak & Tvlendal (12001 1) 
from their population as a whole were £ = 0.27 and e = 0.19. 
The model assumes that entire nebula has the same electron tem- 
perature, T e . The values of T e were calculated based on measure- 
ments of HeII/t4686 A taken fromlTvlenda et al.1 (Tl994l) and the 
method described bv lKalerl dl986l) . With this knowledge we can 
count modeled flux density at 30 GHz. 

Several data sets are required to apply the model de- 
scribed above. He II/ t4686 A line intensities were taken 
Tvlenda et all ([1994), 5 GHz radio flux densities 
Acker etalj (1 19941) . 1.4 GHz flux densities from 
ICond on & Ka plan! (I 19981) and adopted angular diameters 
from ISiodmiak & Tvlendal (120011) . 5 GHz and 30 GHz flux 
density data are from single- dish measurements . This is impor- 
tant because, as noticed by Ziilst ra et al.l d 19891) . observed flux 
densities from a single dish are systematically higher than from 
an interferometer, which is insensitive to the extended, ionised 
halo emission. The data required for the model are available 
for 41 objects out of the 93 for which we have reliable 30 GHz 
flux densities and for these sources the emission model can be 
tested. If source was extended in comparison to OCRA-p beam, 
the adopted angular diameters were used to correct measured 
flux densities. We assumed that PN are spherical sources with 
the given diameter and with constant surface brightness. This 
was a case of only 3 sources: NGC2438, NGC7008, NGC6781. 
Table |2]contains corrected values. 

In Figure [4] we compare the observed and modeled 30 GHz 
flux densities. The model generally agrees well with the mea- 
sured values. The dependence of the logarithm of the observed 
to modeled flux densities at 30 GHz on electron temperature and 
brightness temperature (derived from flux density at 5 GHz - 
S 5 GHz. an d PNe diameter - using T/, — 73.87 • Ssghz/® 2 ) is 
presented in Figures[5]and[6] As can be seen, the differences be- 
tween observed and modeled values are small and independent 
of T e or Tt parameters. 

Most of the data is in good agreement with the model, despite 
its simplicity, so that our observations can be generally explained 
by the free-free emission mechanism with no requirement for an 
additional emission mechanism. However, there are two sources, 
M 2-9 (010.8+18.0) and Vy2-2 (058.3-10.9), in which it can be 
seen that the observed flux density is significantly higher than 
the value expected from the model. 

M 2-9 is a planetary nebula known for its additional emis- 
sion at high er frequencies (see Figu re 7). This emissi on was 
explained by iKwok et al.1 (1 19851) and iDavis et all d 19791) as the 
result of an optically thick stellar wind, which is the dominant 
emission source above about 20 GHz. The high mass loss fr om 
the central star was discovered by Swings & Andrillat (1 19791) on 
the basis of their analysis of its broad Ha line. The observed 
30 GHz point shows that the stellar wind is optically thick above 
30 GHz. 



4 



B. M. Pazderska et al.: Survey of Planetary Nebulae at 30 GHz with OCRA-p 



2.5 



IRAS diagram 



all sources + 
30 GHz measurements o 
Blackbody 



1.5 



0.5 



-0.5 



+ ++ 
-h- + 



"700K 



+ + + + +V+ a V + +++ + v 




100K 



-0.5 



0.5 1 

log 10 (S x (25nm)/S x (12nm)) 



1.5 




Fig. 4. Modeled versus observed flux density at 30 GHz. The 
solid line indicates these flux densities being equal, while the 
dashed line represents the best straight-line fit to relationship be- 
tween log S 3oghz (mod) and log S mem (fibs). 



Fig. 5. Observed to modeled flux density ratio at 30 GHz versus 
electron temperature. 



Vy 2-2 (see Figure 7) is a young and compact nebula. It is 
characterised by bright IR emission, a large optical depth at radio 
frequencies and faint emission. Its turnover frequency is much 
greater than 10 GHz so the assumption of optically thin emission 
at 5 GHz made by our model is not well fulfilled and therefore a 
disagreement with the model is not unexpected. 

Another source which should be mentioned is Hu 1-2 (086.5- 
08.8) (see Figure 7) where emission at 30 GHz is a below the 
expected value. However a check of the other spectral measure- 
ments of this nebula suggest that the point at 5 GHz is overesti- 
mated. 



The model was constructed to check whether any mechanism 
other than free-free emission is needed to expl ain the spectra 
of plan etary nebulae at 30 GHz, as suggested by Casassus et al. 
d2007h . Such an excess was found at 10-30 GHz in the spec- 
tral energy distrib ution (SEP) of sources such as the dark 
cloud LDN 1622 (|Finkbeinerll2004l) . an H II region in Perseu s 
dWatson et alj2005t) and th e Helix nebula dCasassus et al.l2 004). 
iDraine & Lazarianl d 1 998b suggest that radiation from spinning 
electric dipoles created by very small grains (VSGs) or spin- 
ning dust is responsible for the observed excess. Yet in our data 
the effect is not observed. This conclusion is in ag reement with 
new o bservations of planetary nebulae at 43 GHz (lUmana et al.l 
2008), whose measurements are mostly in agreement within er- 
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Fig. 6. Observed to modeled flux density ratio at 30 GHz versus 
brightness temperature. 



rors with those presented in this paper (with exception for PN 
IC2165 (221.3-12.3)). 



6. Conclusions 

We surveyed 442 planetary nebulae from 5-GHz catalogue of 
lAcker et alJ <[l994) with declinations above -15°. 93 sources 
were detected above 23 mJy and their flux densities at 30 GHz 
were measured and presente d in this paper. The mode l for 
free-free emission described by Siodmiak & Tvlendal d2001l) was 
used to predict 30 GHz flux densities based on the 5 GHz obser- 
vations and the known diameters. These calculated values are in 
good agreement with the observed flux densities. No additional 
emission mechanism is needed to explain the observed spectra. 
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Table 2. 30 GHz survey.; col. (1) gives the SESO IAU Planetary Nebula in Galactic coordinates (PNG) name lll.l±bb.b, where 111.1 is Galactic 
longitude and bb.b is Galactic latitude. The values marked by * are sources with optical angular size larger that 1/4 OCRA-p beam, while those 
marked by A are corrected using adopted angu l ar dia meter (see section[5]l.; col. (2) gives the normal descriptive designation, cols. (3) and (4) give 
coordinates for J2000, taken from Acker et al] d!994l) ; cols. (5) and (6) give measured flux densities at 30 GHz and calculated uncertainties. 



PNG 


Name 


OA/ T1AAA\ 

R.A. (J2000) 


Declination (J 2000) 


S 3ocHz (mJy) 


AS 30CHz (mJy) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


a i a o 1 o a 

010.8+18.0 


M 2-9 


17 05 37.952 


1 A AO O A CO 

-10 08 34.58 


63.8 


9.9 


a 1 c n i / Y"> o 

015.9+03.3 


M 1-39 


18 07 30.699 


IO OO A 1 "! tl 1 

-13 28 47.61 


70 


1 1 


A 1 A A AiT O 

019.4-05.3 


M l-6l 


10 A C C C ATI 

18 45 55.072 


1 A cxi—l ci 1—1 AO 

-14 27 37.93 


127 


27 


019.7+03.2 


M 3-25 


18 15 16.967 


1 A 1 A A A A fl 

-10 10 09.47 


76.8 


8.0 


aoa i \ s\ i 1 

020.9-01.1 


M l-5l 


10 O O A/I 1 

18 33 28.943 


1 1 AH O^ A 1 

-11 07 26.41 


229.9 


7.4 


A^O A /"\0 O 

023.9-02.3 


M 1-59 


10 a 1 ^ r\ 1 to 

18 43 20.178 


A A A /I /I O £. A 

-09 04 48.64 


107 


15 


025.3+40.8 


T/— 1 A c 

IC 4593 


1 s 11 A A C A A 

16 ll 44.544 


IO A/1 1^A^ 

+ 12 04 17.06 


66.8 


4.8 


ao co 1 t a 

025.8-17.9 


XT/— 1/— 1 zT O 1 O 

NGC 6818 


1 A /IO n O A A 

19 43 57.844 


1 /I AA 11 A 1 

-14 09 1 1.91 


208 


30 


027.6+04.2 


M 2-43 


18 26 40.048 


AO A O CO y O 

-02 42 57.63 


259 


10 


027.7+00.7 


M 2-45 


1 O in 11 on 

18 39 21.837 


A A 1 A CA AA 

-04 19 50.90 


1 c 

158 


17 


031.0+04.1 


rv 3-6 




1 A A 11 A H 1 A 

+UU 11 4/. IV 


AO A 

98.6 


8.3 


AlO *7 AO A 

U32. /-U2.U 


AA 1 AA 

M 1-00 


1 O CO O/T O/IT 

lo jo zo.z4/ 


A 1 AO /] C -7 A 

-Ul U3 4j. /U 


AO 1 
00. / 


O A 

9.0 


A1/1 A i AO O 

Uj4.U+U2.2 


Is. 3-1 3 


10 AC O/l <CO 


1 r\o r\i 0^ 0^ 
+Uz Ul Id. oj 


OI *2 

91.3 


"7 1 
/.I 


U34.0+1 1.6 


NGC oj 11 


10 n 1.6.% 

lo Iz Uo.3oj 


1 r\/^ ^1 n ni 
+Uu Jl 13. Ul 


OCT 

9o / 


1 Q 
JO 


Al <T 1 AA 1 


A « O 1 

Ap 2-1 


1Q CO 1 f\ A CO 

lo jo 1U.4JV 


1 f\1 0/; C7 1 c 
+U1 30 J /.I J 


O/l 1 

241 


13 


An T *2/l ^ 

U3 /. /-34.J 


XT/" 1 / -1 "7 A AO 

JNuG /UU9 


zl U4 lU.o / / 


11 01 /IC oc. 
-1 1 Zl 4o.zJ 


A *7 A 

4/U 


1 
21 


An C AA Q 

U3 /.o-Uo.3 


NuL /9U 


IV ZZ JO.VDJ 


1 A 1 O A A 11 A C 

+U1 3U 4o.4j 


A*2 

3U3 


1 
12 


A "3 O O i 1 O A 
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Fig. 7. Radio spectra of PN with flux de nsities at 1.4, 5, 14.7 GHz from the literature (see ICondon & Kaplanl 119981 : I Acker et all 1 19941 ; 
Aaquist & Kwok 199Cl lMilne & Alleril 19821) and 30 GHz (OCRA-p) points from Tabled cols. (5) and (6). 



010.8+18.0 (M2-9) 



015.9+03.3 (Ml-39) 



019.4-05.3 (Ml-61) 





+ 







1 10 100 

Frequency (GHz) 



1 10 100 

Frequency (GHz) 



1 10 100 

Frequency (GHz) 



019.7+03.2 (M3-25) 



023.9-02.3 (Ml -59) 



027.6+04.2 (M2-43) 




1 10 100 

Frequency (GHz) 



1 10 100 

Frequency (GHz) 




1 10 100 

Frequency (GHz) 



032.7-02.0 (Ml -66) 



032.7-02.0 (Ml-66) 



034.6+11.8 (NGC6572) 



1 10 100 

Frequency (GHz) 



a 



1 10 100 

Frequency (GHz) 



1 10 100 

Frequency (GHz) 



037.8-06.3 (NGC6790) 



038.2+12.0 (Cn3-1) 



041.8-02.9 (NGC6781) 



1 10 100 

Frequency (GHz) 



1 10 100 

Frequency (GHz) 



1 10 100 

Frequency (GHz) 



045.4-02.7 (Vy2-2) 



046.4-04.1 (NGC6803) 



048.0-02.3 (PB10) 




1 10 100 

Frequency (GHz) 





tr~" * — 


i 




/ 
/ 









1 10 100 

Frequency (GHz) 



100 



1 10 100 

Frequency (GHz) 



B. M. Pazderska et al.: Survey of Planetary Nebulae at 30 GHz with OCRA-p 



Fig. 7. continued. 
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Fig. 7. continued. 
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